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fSJ , Abstract 

In this work we consider a spatially homogeneous and flat Friedmann-Robertson- Walker (FRW) 
space-time filled with non-interacting matter and dark energy components. The equation of state (EoS) 
parameters of the two sources are varied phenomenologically in terms of scale factor of the FRW space- 
" time in such a way that the evolution of the Universe takes place from the early radiation-dominated 

phase to the present dark energy-dominated phase. We constrain the derived model in two cases with 
the latest astronomical observations, and discuss the best fit model parameters in detail. First, we 
qh explore a special case of the model with WMAP+BAO+H0 observations by synchronizing the model 

with the ACDM model at the present epoch. An interesting point that emerges from this observational 
analysis is that the model is not only consistent with the ACDM predictions at the present epoch but 
also is indistinguishable from the ACDM model in revealing the future dynamics of the Universe. In 
the second case, we find observational constraints on general class of the model from Supernova+BAO 
observations. The derived model, in the general case, predicts age of the Universe, Hubble constant, 
density parameters and equation of state parameter of dark energy consistent with the ones obtained 
from seven year WMAP observations. The model advocates cosmological constant as a candidate of 
dark energy, which is consistent with the WMAP observations. Finally, we conclude that the derived 
model offers a unified description of the evolution of Universe from the early radiation-dominated phase 
to the present dark energy-dominated phase in accord with the current astronomical observations. The 
model is physically viable and is applicable to the real Universe. 



Keywords: FRW space-time • Accelerating Universe • Varying equation of state parameters • Dark 
energy ■ Cosmological Constant. 



X: 

1 Introduction 

It is not a matter of debate now whether the Universe is accelerating at the present epoch since it is 
strongly supported by various astronomical probes of complementary nature such as type la supernovae 
data (SN Ia)[U[2], galaxy redshift surveys [3], cosmic microwave background radiation (CMBR) data [HE] 
and large scale structure [6]. Observations also suggest that there had been a transition of the Universe 
from the earlier deceleration phase to the recent acceleration phase [7\. We do not have a fundamental 
understanding of the root cause of the accelerating expansion of the Universe. We label our ignorance 
with the term "Dark Energy" (DE), which is assumed to permeate all of space and increase the rate 
of expansion of the Universe [8J. On the other hand, the inclusion of DE into the prevailing theory of 
cosmology has been enormously successful in resolving numerous puzzles that plagued this field for many 
years. For example, with prior cosmological models, the Universe appeared to be younger than its oldest 
stars. When DE is included in the model, the problem goes away. 



*Department of Applied Mathematics, Delhi Technological University, Delhi-110 042, India. E-Mail: sukuyd@gmail.com 
^Institute of Theoretical Physics, School of Physics & Optoelectronic Technology, Dalian University of Technology, Dalian, 
116024, P. R. China . E-Mail: lxxu@dlut.edu.cn 



1 



The most recent WMAP observations indicate that DE accounts for around three fourth of the total 
mass energy of the universe [9]. However, the nature of DE is still unknown and various cosmological 
probes on theoretical and experimental fronts are in progress to resolve this problem. The simplest 
candidate for the DE is the cosmological constant (A) or vacuum energy since it fits the observational 
data well. During the cosmological evolution, the A-term has the constant energy density and pressure 
Pde = —pde- However, one has the reason to dislike the cosmological constant since it always suffers 
from the theoretical problems such as the "fine-tuning" and "cosmic coincidence" puzzles [10] . That is 
why, the different forms of dynamically changing DE with an effective equation of state (EoS), Wde = 
Pde/ Pde < — 1/3, have been proposed in the literature. Other possible forms of DE include quintessence 
(wde > — 1) [11] ■ phantom (wde < —1) [32] etc. Using the first year Sloan Digital Sky Survey-II (SDSS-II) 
supernova results, Lampeitl et al. [13] found that —0.99 < w^e < —0.65 while the seven year WMAP 
results put Wd e in the range —1.55 < Wde < —0.7 (see, Jarosik et al. [14]). 

In cosmology, the evolution of the Universe is described by the Einstein's field equations together with 
an EoS (p = wp) for the matter content. Usually the field equations are solved and analyzed separately 
for different epochs, i.e., for inflationary phase, radiation-dominated phase and matter-dominated phase. 
Some authors have presented unified solutions for these epochs. For instance, Israelit and Rosen [151 [TB] 
presented a unified EoS, where the pressure varies continuously from pre-matter period (p = —p) to the 
radiation-dominated phase {p = p/3) and then radiation to matter-dominated period (p = 0) for spatially 
closed, flat and open FRW models. They also described a transition between pre-matter to radiation 
and radiation to matter-dominated epoch. Similarly, Carvalho [T7] studied a flat FRW model, where the 
Universe undergoes a transition from an inflationary phase to a radiation-dominated phase. However, 
in the models presented by Israelit and Rosen |15[ I16j . only ordinary matter was taken into account 
while Carvalho [17] studied the model with cosmological constant. In a recent paper [18], a physically 
reasonable and mathematically tractable cosmological model is studied by assuming time-varying EoS 
parameters for matter and DE. Here, we intend to study a cosmological model based on a special case of 
EoS parameters proposed in [15] . 

In this paper, we consider non-interacting matter (dark matter plus baryonic) and DE energy com- 
ponents within the framework of a spatially homogeneous and flat FRW space-time in general relativity. 
Following Ref . |18] , the EoS parameters of the two sources have been varied phenomenologically in terms 
of scale factor of the FRW space-time in such a way that the evolution of the Universe takes place from 
the early radiation-dominated phase to the present DE-dominated phase. The paper is structured as fol- 
lows. The model and field equations are presented in Section [2l In Sectional we present the cosmological 
model with time-varying EoS parameters for the matter and DE components. In Section U we explore 
a special case of the model with WMAP+BAO+H0 observations by synchronizing the model with the 
ACDM model at the present epoch. Section [5] is devoted to find the best fit model constrained with latest 
SN Ia+BAO observations. The findings of the paper are summarized in Section [HI 



2 Model and field equations 

We consider a spatially homogeneous and flat FRW line element that applies to the real Universe. It is 
written as 

ds 2 = -dt 2 + a 2 (t) [dr 2 + r 2 (dd 2 + sin 2 ## 2 )] , (1) 

where a(t) is cosmic scale factor and other symbols have their usual meanings. 

The Einstein's field equations in case of a mixture of matter and DE components, in the units 8nG = 
c = 1, read as 

Rfiv — ~^9fivR = ~Tfivi (2) 
where T^ u = Tffl + T^f 1 is the overall energy momentum tensor with Tjffl and Tjffl as the energy 
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momentum tensors of matter and DE, respectively. These are given by 

T (m) n _ diag [_ /£ , TO) Pm iPmi Pm ] 

= diag [-1, w m ,w m , w m ]p m (3) 

and 

T (de) n = diag [_ pdfti pd£ ^ pde: p<fe ] 

= diag [-1, W de ,W de , W de ]Pde (4) 

where p m and p m are, respectively the energy density and pressure of the matter fluid while w m = p m /p m 
is its EoS parameter. Similarly, pd e and pd e are, respectively the energy density and pressure of the DE 
fluid while Wde = Pde/Pde is the corresponding EoS parameter. 

In a comoving coordinate system, the field equations ([2]) for the space-time ([1]) , in case of © and (|3J) , 
read as 

2H + 3H 2 = -W m p m - WdePde, (5) 

3H 2 = p m + p de - (6) 

Here an over dot indicates ordinary derivative with respect to t, and H = a/ a is the Hubble parameter. 
The energy conservation equation T ;„ = yields 

p m + 3(1 + w m )p m H + p de + 3(1 + w de )pdeH = 0. (7) 

We assume that the matter and DE components are non-interacting. Therefore, the energy momentum 
tensors of the two sources may be conserved separately. 

The energy conservation equation T ^ v . v = 0, of the matter fluid leads to 

p m + 3(1 + w m )p m H = 0, (8) 
whereas the energy conservation equation T^ de ^ ^ v - v = 0, of the DE component yields 

Pde + 3(1 + W de )pdeH = 0. (9) 

3 Cosmology with varying EoS parameters 

In order to construct a model for unified description of the evolution of the Universe from the early 
radiation-dominated phase to the recent DE-dominated phase, we assume the following forms for the EoS 
parameters of matter and DE (see [18] for more general EoS parameters): 

Wm = J^hry (10) 

wx a 

Wde = a . : , 11 
x a + 1 

where x = a /a* with a* being some reference value of a. Further, a is some positive constant parameter 
while w is a negative constant. 

Substituting ([TO]) into ©, we find 

Pm = ~ 4 , (12) 
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~ 3x 4 ' { 1 



where C\ is a positive constant of integration. 
Similarly, substituting (fTT|) into ([9]), we obtain 



a ~ ° 2 ^ + 1) ~ W/a ri4l 

= ^) ' ( 15 ) 

where C 2 is a positive constant of integration. 
• For x <C 1, we have 

^ Ci _ 1 _ 

4 j Pm — „ Pm i Pde ~ o j 
6 X 3 

which means that for x <C 1 matter dominates over DE, as expected. Besides, in this limit we find 



wC 2 

Pde 



3x( 3 - Q ) ' 
For x S> 1, we get 



3x (3+a) • 



Since a > 0, the matter pressure decreases with the expansion of the Universe much faster than the 
density. Thus, p m » p m as required by a matter dominated Universe. Note that p m decreases exactly as 
in cosmological dust models. 

In addition, in this limit we obtain 

C 2 

Pde=W Pde: p de ^-^—^. 

Since w < 0, it follows that pd e S> p m and \pde\ S> p m for s > 1. 
Now from (El), it follows that 



H _ / Cl( ^ + i)'/° + C2 (x» + ir^ 

V 3ar 3x d 

Using (JT0]) - (fT2|) . ([T4"j) and (fT6j) we find that © is satisfied identically, as expected. For a detailed classical 
treatment of the above equations in case of more general EoS parameters, the reader is referred to Ref. 
|18j . Here, we are interested in finding the observational constraints on the cosmological model in hand. 
The effective EoS parameter is obtained as 

Peff Pm+Pde 3C 2 WX 1+a + C\ (1 + X a ) 

W e ff = = = r 1 + 3-1 • 

Peff Pm+Pde 3(1 + X a ) C 2 X + Cj ( 1 + X a ) a 

In terms of the cosmological redshift z we can write a = ao/(l + z), where ao is the present day value 
of a, which corresponds to z = 0. Thus, 

1 + Z * fia\ 

X = TT7- (18) 

It follows that x varies from to oo as z varies from oo to —1. 

In the following section, we consider a special class of the proposed model and compare it with the 
ACDM model by subjecting the model to WMAP+BAO+H0 observations. 
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4 Observational constraints on a special class of the model 

Table [T] shows the variation of the EoS parameters as z varies from oo to — 1 . 



Table 1: Extreme values of EoS parameters 



EoS Parameter 


z — > oo 


z-t-1 


Wm 


1/3 










;7; 


W eff 


1/3 


w 



It is interesting to observe that the constant w is decisive in future dynamics of the Universe. Also 
this constant is at our discretion. Thus, if we choose —1 < w < —1/3, the dynamically evolving DE will 
never cross the phantom divide line (PDL) {wde = — 1) and the accelerated expansion of the Universe 
will continue in future with quintessence form of DE. If we set w < —1, the PDL will be crossed and 
the Universe will enter the phantom regime. Choosing w = —1, we ensure that dynamics of the future 
Universe will be purely governed by cosmological constant. In what follows we shall carry on with the 
choice w = — 1. Also, in the proposed model it seems reasonable to interpret a* as the value of a for 
which p m = p de - Then, from ([12]) and ([33} it follows that C x = C 2 2 2/a . 

Taking into account the above considerations and restoring the SI units with Mp = hc/8-KG (Mp 
being the reduced Planck mass), equations (fT2]) - (fT7|) reduce to 



C 2 2 2 /° (x a + 1) 



1/q 



(19) 



Pr, 



C 2C 2 2 2/a (^+1)0-°OA» 

3x 4 

_ C 2 (x a + lf a 

Pde — q j 



2 /V" _i_ -n(3-a)A* 



Pde 



C 2 c 2 {x a + 1 



,(3-«) 



(20) 
(21) 
(22) 



H 



I c 2 hc 

3Mpx 4 



2 2/a ( x a + !)!/« + x + X fh 

3x 1+a + 2 2 / a (x a + l)- 2 /° 



UK 



eff 3(x« + l) [x + 2 2 /«(x« + l)-2/«]' 
The deceleration, jerk and snap parameters of the model are respectively given by 

a 



q = - 



. H' _ 2 2 / a {x a + 2) - x{2x a - \){x a + 1) 2 / Q 
aH 2 ~ ~ X ~H~ 2{x a + 1) [2 2 /« + x(x a + l) 2 / a ] 



a H n 

—3=1 + 4x— + x 2 
an 6 H 




(23) 
(24) 

(25) 
(26) 
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s = 



aH 4 



a 




(27) 




(28) 



3Mj,H 2 2 2 A* + x(x a + l) 2 /« ' 



x(x a + l) 2 /' 



(29) 



de 



3M 2 # 2 



2 2/q + + • 



In view of ([6]), it is observed that Q m + f^e = 1- 

In order to examine where the unified model stands with respect to the so called standard ACDM 
model, we produce the kinematics of the standard ACDM model in Appendix I. 

In what follows, we find observational constraints on the parameters of the derived model and ACDM 
model using the following observational results from WMAP7+BAO+-£fo given in Ref.|14j: 



the constraints on deceleration parameter of the ACDM model are q^o = — 0.59_ 02 . WMAP observations 
suggest that the flat ACDM model is most suitable to describe the Universe at the present epoch. So 
it would be useful to synchronize the derived model with the values of Hubble parameter, deceleration 
parameter and DE density parameter as given above. With the best fit values of dark energy density 
parameter and deceleration parameter, we immediately find a = 18.15 and z* = 0.42. It may be noted 
that the values of a and are highly sensitive with respect to the values of Qde and q in their permissible 
domains from WMAP. For instance, choosing QdeO = 0.743 and c/ao = —0.57, we obtain a = 8.06 and 
z* = 0.50. However, in what follows, we shall utilize only the best fit values of the parameters derived 
above. 

Table [2] displays the extreme values of various parameters of the derived model along with their best 
fit values in contrast with those of the ACDM model. For the sake of completeness and to examine how 
the derived model differs from the ACDM model, we show the variation of various parameters for the 
derived model and the ACDM model in Fig. [I] to Fig. [8l After a careful and straightforward analysis of 
the values of various parameters displayed in Table [2] and the graphics (Fig. Q] to Fig. [8|) , we conclude 
that the model in hand is not only consistent with the ACDM model at the present epoch but also is 
indistinguishable from the ACDM model in revealing the future dynamics of the Universe. Next, it differs 
significantly from the ACDM model at the earlier epochs of evolution dictating its advantages over the 
usual ACDM one. For instance, the ACDM model accounts only for pressureless matter even in the early 
stages of evolution. On the other hand, the derived model successfully describes the evolution of the 
Universe from the early radiation-dominated matter phase to the current pressureless matter phase (see 
Fig. [T] and values of w m in Table [2]) . One may also observe the significant difference in the values of the 
parameters q, j and s for the two models (see Table [2] and Fig. [6] to Fig. [8]). 




m 



2 while 
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Table 2: Extreme and the best fit values of various parameters pertaining to the derived model and the 
ACDM model 



Model ->■ 


Derived Model 


ACDM 


Parameter 


z — y oo 


z = 


z -> -1 


z — > oo 


z = 


z — >- — 1 


Wrr, 


1/3 


0.0005 



















-0.998 


-1 


-1 


-1 


-1 


Wf, f f 


1/3 


-0.726 


-1 





-0.728 


-1 


Vrr, flCT 10 Pa) 


OO 


0.001 














VriP flO~ 1(J Pa) 





-6.073 


-6.081 


-6.083 


-6.083 


-6.083 


Vpff flO~ lu Pa) 


oo 


-6.072 


-6.081 


-6.083 


-6.083 


-6.083 


Pm (lO -27 kg m -3 ) 


oo 


2.525 





oo 


2.525 





Pri P (10~ 27 kg m -3 ) 


oo 


6.759 


6.757 


6.759 


6.759 


6.759 


p P ff (10~ 27 kg m -3 ) 


oo 


9.284 


6.757 


oo 


9.284 


6.759 




1 


0.272 





1 


0.272 










0.728 


1 





0.728 


1 




1 


1 


1 


1 


1 


1 


# (km s- 1 Mpc" 1 ) 


oo 


70.4 


60.1 


oo 


70.4 


60.1 


Q 


1 


-0.59 


-1 


0.5 


-0.59 


-1 


3 


3 


1.03 


1 


1 


1 


1 


s 


-15 


-0.79 


1 


-3.5 


-0.22 


1 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

Z 



Figure 1: EoS parameters w m (dashed blue curve), WAm (solid blue curve), Wde (dashed red curve), w^de (solid red curve) 
, w e ff (dashed green curve), u>Ae// (solid green curve) vs z. The vertical dashed line is for z — 0. 
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Figure 2: Pressures p m (dashed blue curve), PAm (solid blue curve), pde (dashed red curve), pAde (solid red curve) , p R ff 
(dashed green curve), PAeff (solid green curve) vs z. The vertical dashed line is for z — 0. The curves related to PAde and 
PAeff are overlapping as p A de = PAeff for all z. 
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Figure 3: Energy densities p m (dashed blue curve), pAm (solid blue curve), pde (dashed red curve), pAde (solid red curve) , 
p e ff (dashed green curve), PAeff (solid green curve) vs z. The vertical dashed line stands for z = 0. 




Figure 4: Density parameters Q m (dashed blue curve), fiAm (solid blue curve), ttde (dashed red curve), S^Ade (solid red curve) 
, 0.eff (dashed green curve), Q,Aeff (solid green curve) vs z. The vertical dashed line is for z = 0. Since f2 e // = fW// = 1 
for all z, the green curves related to these parameters are overlapping. 



8 




Figure 5: Hubble parameters H (dashed blue curve), Ha (solid blue curve) vs z. The vertical dashed line is for 2 = 0. 
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Figure 6: Deceleration parameters q (dashed blue curve), qa (solid blue curve) vs z. The horizontal and vertical dashed 
lines respectively stand for q — and z = 0. The transition redshift for the derived model is zt = 0.35 while for the ACDM 
model the transition takes place at zt = 0.75. 
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Figure 7: Jerk parameters j (dashed blue curve), jh (solid blue curve) vs z. The vertical dashed line stands for z = 0. The 
jerk parameter j of the derived model attains its maximum value 5.42 at z = 0.43. 
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Figure 8: Snap parameters s (dashed blue curve), sa (solid blue curve) vs z. The horizontal and vertical dashed lines 
respectively stand for s — and z = 0. The snap parameter s of the derived model attains its maximum value 7.67 at 
2 = 0.54 and minimum value —37.26 at z = 0.33. 



5 Observational constraints on general class of the model 

The Friedmann equation of the derived model, in the general case, can be rewritten as 

8ttG 



{(q° + q?)/(l + a ?)} 1 / a {(q" + <)/(l + a»)}- 3 "/ a 

PmO 7 r PdeO q 



Hi 



m0" 



{(a a + a«)/(l + a?)} 



l/a 



+ 0, 



deO" 



{(a^ + a^/a + a-)}" 3 ^ 



where 







deO 



1. 



(30) 



(31) 



and a is the scale factor. 

To test the viability and to obtain the parameter space of this model, we use SN la and BAO data 
sets and the Markov Chain Monte Carlo (MCMC) method. Our code is based on the publicly available 
package cosmoMC [24J. At first, we modified the code to add three new parameters a, w and a*. The 
following 4-dimensional parameter space is adopted 



P = {w c ,a,w,a*} 



(32) 



where w c = Q c h 2 is the physical cold dark matter density. We take the following priors to model 
parameters: u c G [0.01,0.99], a E (0, 100], a* G [0,0.1] and w £ [—3,0]. In addition, the hard coded prior 
on the comic age lOGyr < to < 20Gyr is imposed. Also, we fixed the physical baryon density a>b = 0.022 
|25j from big bang nucleosynthesis and the new Hubble constant H = 74.2 ± 3.6 km s" 1 Mpc" 1 [26] . 



To get the distribution of parameters, we calculate the total likelihood C oc e x / 2 , where \ z is given 



as 



X 2 = Xbao + Xsn- (33) 

The 557 Union2 data [27] with systematic errors and BAO [281129] are used to constrain the background 
evolution. For the detailed description, see Refs. [30] . 

After running 8 independent chains and checking the convergence to stop sampling when the worst 
e-values [the variance(mean)/mean(variance) of 1/2 chains] R — 1 is of the order 0.01, the global fitting 
results are summarized in Table [3] and Fig. [9] For the sake of comparison and to see the viability 
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of the derived results, we also show the values of various parameters in Table [3] based on WMAP and 
WMAP+BAO+HO observations (see Jarosik et al. [J3]). 

One may see that the derived model is in close agreement with the results predicted by WMAP and 
WMAP+BAO+HO observations. It may be noted that the values of a can be taken in a large range, and 
current cosmic observations from SN la and BAO cannot give a tight constraint to the parameter a. For, 
smaller values of a* and values of a > 1, in view of (I30p . lead to the DE model with EoS w = w via 

H 2 (a) « Hi [n m0 a- 3 + n^a^ 1 ^] • (34) 

Further, in view of (1101) it deserves mention that smaller values of a* and values of a > 1 yield values of 
w m closer to 0. This is turn implies that the observations from SN la and BAO force the model to describe 
the evolution of the Universe from relatively later phase of radiation to the current phase dominated by 
some sort of DE with EoS w ~ — 1. Thus, the model puts forward cosmological constant as a candidate 
of DE. This is consistent with the WMAP observations. 



Table 3: Best fit values of the derived model parameters along with error bars at 68% and 95% levels, 
and values of some relevant parameters based on WMAP and WMAP+BAO+HO observations 



Parameter 


SN Ia+BAO 
(This paper) 


WMAP 

(Jarosik et al. [Til) 


WMAP+BAO+HO 
(Jarosik et al. [13]) 


a 


A 9 009+15. 164+37.349 
* z - ooz -42.832-42.832 






a* 


n nnni qp+0.000038+0.000109 

U.UUU±OO_ .000138-0 .000138 








fl 9Q7+0.0200+0.0427 
U.zo < _ .0197-0 .0368 






n c h 2 


n 1 qo+0.0115+0.0254 
U.±dO_ .0112-0.0188 


u.±iuy_ 00056 


n 1 19Q+0.0035 

u.±izo_Q.oo 35 


Q A 


n 71 q+0.0197+0.0368 
u - ' lo -0.0200-0.0427 


704+0.029 
u -'°*-0.029 




w 


1 n7^s+°- 0943 + - 186 

J..U/OO_ .0944-0 .186 


1 1 9 +0.42 
1 - LZ -0.43 


-0.9801H1 


Age (Gyr) 


1011 q+0. 499+0.640 
lO.llJ_ 4so_o .983 


13.75i°il 


13.75±g:g 


H Q (km s _1 Mpc" 1 ) 


70 fi1 7+2.870+5.990 
'°- ui ' -2.775-5.0607 


n.o±3:| 
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Figure 9: The ID marginalized distribution on individual parameters and 2D contours with 68% and 
95% confidence limits obtained by using SN Ia+BAO data points. The shaded regions show the mean 
likelihood of the samples. 

6 Concluding remarks 

In this work, we have investigated a cosmological model within the framework of a spatially homogeneous 
and flat FRW space-time filled with non-interacting matter and DE components. Following Ref. |18j . the 
model is derived by assuming time- varying EoS parameters of the two sources, which in turn provides an 
elegant evolution of the Universe from the early radiation-dominated phase to the present DE-dominated 
phase. We have explored a special case of the model with WMAP+BAO+HO observations by synchroniz- 
ing the model with the ACDM model at the present epoch. The observational analysis suggests that the 
model is not only consistent with the ACDM predictions at the present epoch but also is indistinguishable 
from the ACDM model in revealing the future dynamics of the Universe. Thus the derived model, in the 
special case, has already revealed what ACDM model has to offer. In addition, it accounts for radiation- 
dominated matter phase at early epochs. Thus, the derived model has an advantage over the usual 
ACDM one. We have also tested the viability of the general class of the model by constraining the model 
with SN la and BAO data sets. In the general case also the derived model yields parameters consistent 
with the WMAP and WMAP+BAO+HO observations. The model advocates cosmological constant as a 
candidate of DE, which is consistent with the WMAP observations. Finally, we conclude that the derived 
model offers a unified description of the evolution of Universe from the early radiation-dominated phase 
to the present DE-dominated phase in accord with the current astronomical observations. The model is 
applicable to the real Universe, and is supposed to yield more accurate results with the advancement of 
cosmic data. It would be interesting to find observational constraints on the cosmological model based 
on generalized EoS parameters for matter and DE reported in Ref. [18| . 
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Appendix I. Elements of ACDM cosmology 

The standard ACDM Universe is governed by the scale factor 



2 



a\ = ai sinhi l^J—^-tj , (35) 
where a\ is a constant. 

Using the relation a a = ao/(l + z), the Hubble parameter, deceleration parameter, jerk and snap 
parameters in terms of the redshift are obtained as 



H A = ^ W( y + a8] , (36) 

-2fl3 + q3(l + z)3 

QA ~ 2[o§ + a?(l + z) 3 ]' [67) 

3 A = 1, (38) 



2a 3 1 -7a?(l + z) 3 
2[ag + af(l + 2i) 

The pressure and energy density of the ordinary matter in ACDM cosmology are 



SA = 70 ?; (39) 



PAm = 0, (40) 

MlKca\(l + zf 

PAm = jfcg ' ( 41 > 

while the pressure and density of the vacuum energy associated with the cosmological constant read as 

M P A(3 

PAde = ^— , (42) 

M 2 p Kc 

PAde = —J—- (43) 

The density parameters and effective EoS parameter in the ACDM cosmology is given by 

qf(l + z) 3 , , 

UAm ~ a 3 + a3(l + z)3]' ^ 
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""■- Wi+^r (45> 

w --" = -^TWTWY <46) 

Obviously, the EoS parameters of ordinary matter and vacuum energy in the ACDM Universe are 
WAm = and WAde = — 1) respectively. 
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